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Two classes of gamma-ray bursts were identified in the BATSE catalogs characterized by their 
durations. There were also some indications for the existence of a third type of gamma-ray bursts. 
Swift satellite detectors have different spectral sensitivity than pre-Swift ones for GRBs. Therefore 
in this paper we analyze the bursts' duration distribution and also the duration-hardness bivariate 
distribution, published in The First BAT Catalog. Similarly to the BATSE data, to explain the BAT 
GRBs' duration distribution three components are needed. Although, the relative frequencies of the 
groups are different than they were in the BATSE GRB sample, the difference in the instrument 
spectral sensitivities can explain this bias in a natural way. This means theoretical models may have 
to explain three different type of gamma-ray bursts. 



I. INTRODUCTION 

The discovery of the third type of GRBs goes back 
as early as 1998 [l|, After that many research 
groups studied the BATSE bursts' sample and con- 
cluded the third group of the GRBs statistically ex- 
ists [1, 0, [E H, 0, S H, ■ Later several pap ers were 
published analysing different data sets [ill Il2j . In 
this paper we analyze the bursts' duration distribu- 
tion and also the duration-hardness bivariate distri- 
bution, published in The First BAT and The BATSE 
Catalog. 



N 



L = log? = J^log (g(xi,pi, ...,p k )) 



(2) 



The ML procedure maximizes L according to 
pi, ■■■,Pk- Since the logarithmic function is monotonic 
the logarithm reaches the maximum where / does it as 
well. The confidence region of the estimated param- 
eters is given by the following formula, where L max 
is the maximum value of the likelihood function and 
I/o is the likelihood function at the true value of the 
parameters [HI]: 



II. THE ONE DIMENSIONAL GAUSSIAN 
FITS 

We used the Maximum Likelihood (ML) method for 
the analysis of the (Swift) BAT and BATSE bursts. 
The ML method assumes that the probability den- 
sity function of an x observable variable is given in 
the form of g{x,p\, ■ ■■,Pk) where p\, ...,pk are param- 
eters of unknown value. Having N observations on x 
one can define the likelihood function in the following 
form: 



JY 



1 = Y\9{x t ,pi, ...,Pk) 



(1) 



or in logarithmic form (the logarithmic form is more 
convenient for calculations): 



2(L r 



L ) 



xt 



(3) 



Therefore one can fit the logTgo distribution using 
ML with a superposition of k log-normal components, 
each of them having 3 unknown parameters to be fit- 
ted with N measured points. Our goal is to find the 
minimum value of k suitable to fit the observed dis- 
tribution. Assuming a weighted superposition of k 
log-normal distributions one has to maximize the fol- 
lowing likelihood function: 



N 



L k = ^ log ^wiftixi, log Ti,a{) 



(4) 



where wi is a weight, fi a log-normal function with 
log Ti mean and c; standard deviation having the form 
of 
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FIG. 1: The duration distribution of the BATSE bursts and the two- and three-Gaussian fits. 



, 1 / (x — logTj) 2 
J I = 7= ex P o~2 

and due to a normalization condition 



J2 w i = N 



(5) 



(6) 



1=1 



with three functions with eight parameters (three 
means, three standard deviations and two weights). 
The logarithm of the best likelihood (L3) is 12326.25. 
According to the mathematical theory, twice the dif- 
ference of these numbers follows the x 2 distribution 
with three degrees of freedom because the new fit has 
three more parameters. The difference is 6.14 which 
gives us a 0.5% probability. Therefore the three- 
Gaussian fit is better and there is only a 0.005 prob- 
ability that it is caused by statistical fluctuation. 



III. ONE DIMENSIONAL GAUSSIAN FITS 
FOR THE BATSE BURSTS 

A fit to the duration distribution of the BATSE 
bursts was taken using a maximum likelihood method 
with the superposition of two log-normal distribu- 
tions. This can be done by a standard search for 5 
parameters with N = 1929 measured points. Both 
log-normal distributions have two parameters; the 
fifth parameter defines the weight (101) of the first 
log-normal distribution. The second weight is w 2 = 
(N — w\) due to the normalization. Therefore we ob- 
tain the best fit to the 5 parameters through a maxi- 
mum likelihood estimation. The maximum of the like- 
lihood was 7,2=12320.11. 

Secondly, a three-Gaussian fit was made (Figure 1. 
shows both the two and the three component fits) 



IV. ONE DIMENSIONAL GAUSSIAN FITS 
FOR SWIFT GRBS 

In the Swift BAT Catalog [3 there are 237 GRBs, 
of which 222 have duration information. Fig. 2. shows 
the logTgo distribution. We made fits for this distri- 
bution. Assuming only one log-normal component the 
fit gives Li max — 951.666 but in the case of k=2 one 
gets L 2m ax = 983.317 . 

Based on Eq. §5§ we can infer whether the addi- 
tion of a further log-normal component is necessary 
to signifincantly improve the fit. We make the null 
hypothesis that we have reached already the the true 
value of k. Adding a new component, i.e. moving from 
k to k + 1, the ML solution of L kmax has changed to 
L(k+i)maxi but Lq remained the same. In the mean- 
time we increased the number of parameters with 3 
(tOfe+i, logT k+ i and <7( fe+ i)). Applying Eq. © on 
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FIG. 2: The duration distribution of the Swift bursts and a three-Gaussian fit. 



both Lkmax and L/^+^max we § e t after subtraction 

2{L(k+l)max ~ Lkmax) ~ x\ (7) 

For k — 1 L2rnax is greater than L\ max by more than 
30, which gives for x| an extremely low probability of 
5.88 x 1CP 13 . It means the two log-normal fit is really 
a better approximation for the duration distribution 
of GRBs than one log-normal. 

Thirdly, a three-log-normal fit was made combin- 
ing three /& functions with eight parameters (three 
means, three standard deviations and two weights). 
The highest value of the logarithm of the likelihood 
{L^ max ) is 989.822. For two log-normal functions the 
maximum was L^max — 983.317. The maximum thus 
improved by 6.505. Twice of this is 13.01 which gives 
us the probability of 0.46% for the difference between 
Limax and Lj, max is being only by chance. Therefore 
there is only a small chance the third log-normal is not 
needed. Differently said the three-log-normal fit (see 
Figure 2.) is better and there is a 0.0046 probability 
that it was caused only by statistical fluctuation. 

One should also calculate the likelihood for four 
log-normal functions. The best logarithm of the ML 
is 990.323. It is bigger with 0.501 than it was with 
three log-normal functions. This gives us a low signif- 
icance (80.1%), therefore the fourth component is not 
needed. 



V. THE TWO DIMENSIONAL GAUSSIAN 
FITS 

When studying a GRB distribution, one can assume 
that the observed probability distribution in the pa- 
rameter space is a superposition of the distributions 
characterizing the different types of bursts present in 
the sample. Using the notations x and y for the vari- 
ables (in a 2D space) , and using the law of full prob- 
abilities, one can write 

fc 

P(x,y) = ^2p(x,y\l)pi (8) 
l=i 

In this equation p(x, y\l) is the conditional probability 

density assuming that a burst belongs to the i-th class. 

Pi is the probability for this class in the observed sam- 
k 

pie Pi — 1); where k is the number of classes. In 
1=1 

order to decompose the observed probability distribu- 
tion p(x, y) into the superposition of different classes 
we need the functional form of p(x,y\l). The proba- 
bility distribution of the logarithm of durations can 
be well fitted by Gaussian distributions, if we restrict 
ourselves to the short and long GRBs [2j . We assume 
the same also for the y coordinate. With this assump- 
tion we obtain, for a certain i-th class of GRBs, 
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FIG. 3: The duration - hardness distribution of the BATSE bursts and the three components. 



P(x,y\l) 



2iro~ x o~ y y/l — r 2 



x exp 



2(1 



(x - a x ) 2 (y-a y ) 2 2r(x - a x )(y - a y ) 



' x u y 



(9) 



where the means, <r x , o~ y are the dis- 

persions, and r is the correlation coefficient. Hence, 
a certain class is defined by 5 independent parame- 
ters, a x , a y , <J X , o~ v , r, which are different for different 
I, If we have k classes, then we have (6k — 1) inde- 
pendent parameters (constants), because any class is 
given by the five parameters of Eq. ^ and the weight 

pi of the class. One weight is not independent, be- 

fe 

cause ^2 Pi = 1 • The sum of k functions defined by 

i=i 

Eq.© gives the theoretical function of the fit. 

By decomposing p(x, y) into the superposition of 
p(x, y\l) conditional probabilities one divides the orig- 
inal population of GRBs into k groups. Decomposing 
the left-hand side of Eq. ([5]) into the sum of the right- 
hand side, one needs the functional form of p(x,y\l) 
distributions, and also k has to be fixed. Because we 
assume that the functional form is a bivariate Gaus- 
sian distribution (see Eq.©), our task is reduced to 
evaluating its parameters, k and pi. 

Balazs et al. [l5| used this method for k = 2, and 
gave a more detailed description of the procedure. 
However, that paper used fluence instead of hardness, 
and used the BATSE data. Here we will make similar 



calculations for k — 2, k 
and Swift observations. 



3 and k = 4 using CGRO 



VI. THE TWO DIMENSIONAL GAUSSIAN 
FITS FOR BATSE BURSTS 

The data from The Final BATSE GRB Catalog 
have been used, in which there are 2702 GRBs, for 
1956 of which both the hardnesses and durations are 
measured. These 1956 GRBs define the sample stud- 
ied in here. The hardness duration distribution can 
be seen in Figure 3. 

Moving from k = 2 to k = 3 the number of pa- 
rameters m increases by 6 (from 11 to 17), and L max 
grows from 1193 to 1237. Since X17 = Xii + X§ the 
increase in L max by a value of 44 corresponds to a 
value of 88 for a Xe distribution. The probability for 
xi — 8^ i s extremely low (< 10 -10 ). Therefore we 
may conclude that the inclusion of a third class into 
the fitting procedure is well justified by a very high 
level of significance. 

Moving from k — 3 to k ~ 4, however, the improve- 
ment in L max is only 6 (from 1137 to 1143) corre- 
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FIG. 4: The duration - hardness distribution of the Swift bursts and the three components. 



sponding to \% > 12, which can happen by chance 
with a probability of 6.2 %. Hence, the inclusion of 
the fourth class is not justified. We may conclude from 
this analysis that the superposition of three Gaussian 
bivariate distributions - and only these three ones - 
can describe the observed distribution of the BATSE 
data Figure 3. also shows the three components. 

VII. THE TWO DIMENSIONAL GAUSSIAN 
FITS FOR THE SWIFT BURSTS 

In the Swift BAT Catalog [3 there are 237 GRBs, 
of which 222 have duration information. Following the 
same procedure of data reduction we have extended 
this sample with all the bursts detected until mid De- 
cember 2008 (ending with GRB 081211a). Our total 
sample thus comprises the first four years of the Swift 
satellite (since the detection of its first burst GRB 
041217) and includes 342 bursts. 222 from Sakamoto 
et al. [14| and 120 reduced by us. The data reduction 



was done by using HEAsoft v. 6. 3. 2 and calibration 
database v. 20070924. For lightcurves and spectra we 
ran the batgrbproduct pipeline. We fitted the spec- 
tra integrated for the duration of the burst with a 
power law model and a power law model with an ex- 
ponential cutoff. As in Sakamoto et al. we have 
chosen the cutoff power law model if the \ 2 of the fit 
improved by more than 6. 

For calculating the hardness ratio we have chosen 
fluence 2 (25-50fceF) and fluence 3 (50-100/ceV) and 
the hardness is defined by the HR = F3/F2 ratio. 
The duration - hardness distribution can be seen in 
Figure 4. The two dimensional fits have been made in 
this plane. 

Moving from k — 2 to k — 3 the number of parame- 
ters m increases by 6 (from 11 to 17), and L max grows 
from 506.6 to 531.4. Since X17 = Xii+Xe the increase 
in L max by a value of 25 corresponds to a value of 50 
for a x\ distribution. The probability for \\ > 50 is 
very low (10~ 8 ), therefore we conclude that the inclu- 
sion of a third class into the fitting procedure is well 
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justified by a very high level of significance. 

Moving from k = 3 to k — 4, however, the improve- 
ment in L max is 3.4 (from 531.4 to 534.8) correspond- 
ing to Xg > 6.8, which can happen by chance with 
a probability of 33.9 %. Hence, the inclusion of the 
fourth class is not justified. We conclude again that 
the superposition of three Gaussian bivariate distri- 
butions - and only these three ones - can describe the 
observed distribution of the Swift BAT bursts. 

One can see from the fits and olso from Figure 4. 
that the mean hardness of the intermediate class is 
very low - the third class is the softest one. This is 
in a good agreement with the BATSE fits (for more 



details see Horvath et al. [9j), where we found that 
the intermediate duration class is the softest in the 
BATSE database. In that database 11% of all GRBs 
belonged to this group. In our analysis p2 — 0.296, 
therefore 30% of the Swift bursts belong to the third 
(intermediate duration) group. 
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